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We study the statistical properties of the scattering matrix S(q\k) for the problem of the scattering 
of light from a randomly rough one-dimensional surface, defined by the equation X3 = £(xi), where 
the surface profile function ((xi) constitutes a zero- mean, stationary, Gaussian random process, 
through the effects of S(q\k) on the angular intensity correlation function C(q, k\q' , k'). The existence 
^—K ' of both the C ^ and C^ 10 - 1 correlation functions is consistent with the amplitude of the scattered 

, field obeying complex Gaussian statistics in the limit of a long surface. We show that the deviation 

of the statistics of the scattering matrix from circular Gaussian statistics and the C^ 10 ' correlation 
function are determined by exactly the same statistical moment. As the random surface becomes 
rougher, the amplitude of the scattered field no longer obeys complex Gaussian statistics but obeys 
complex circular Gaussian statistics instead. In this case, the C*- 10 ' correlation function should 
vanish. This result is confirmed by numerical simulation calculations. 
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l.INTRODUCTION 

The scattering of light from randomly rough surfaces has attracted attention over many years. The majority of 
the theoretical and experimental studies of such scattering has been devoted to coherent interference effects occuring 
in the multiple scattering of electromagnetic waves from randomly rough surfaces and the related backscattering 
O , enhancement phenomenon. 

Recently, attention has begun to be directed toward theoretical [1-12] and experimental [2,7,8,12,13] studies of 
multiple-scattering effects on higher moments of the scattered field, in particular on angular intensity correlation 
functions. These correlation functions describe how the speckle pattern, formed through the interference of randomly 
scattered waves, changes when one or more parameters of the scattering system are varied. 

The interest in these correlations has been stimulated by the expectation that, just as the inclusion of multiple- 
scattering processes in the calculation of the angular dependence of the intensity of the light that has been scattered 
_ ■ incoherently from, or incoherently through, a randomly rough surface, led to the prediction of enhanced backscattering 
ON ! [14] and enhanced transmission [15], their inclusion in the calculation of higher-order moments of the scattered or 
transmitted field would also lead to the prediction of new physical effects. This expectation was prompted by the 
results of earlier theoretical [16,17] and experimental [18-20] investigations of angular intensity correlation functions 
in the scattering of classical waves from volume disordered media. In a theoretical investigation [9] it was predicted 
that three types of correlations occur in such scattering, viz. short-range correlations, long-range correlations, and 
*^ infinite-range correlations. These were termed the C^ x \ C^ 2 \ and correlations, respectively. The correlation 
function includes both the "memory effect" and the "reciprocal memory effect" [9,10], so named because of the 
wave vector conservation conditions they satisfy. Both of these effects have now been observed in volume scattering 
experiments [16,17]. The C^ 2 ' correlation function has also been observed in volume scattering experiments [18,19], 
as has the correlation function [20]. 

Until recently, only the correlation function had been studied theoretically and experimentally [1-8]. In a recent 
series of papers devoted to theoretical studies of angular correlation functions of the intensity of light scattered from 
one-dimensional [9,10] and two-dimensional [10] randomly rough metal surfaces the long-range and infinite 
range correlation functions were calculated, and two additional types of correlation functions, a short-range 
correlation function, named C^ 10 \ and a long-range correlation function, named C^' 5 \ correlation functions were 
predicted. In very recent experimental work [12] the envelopes of the and C^ 10 ^ correlation functions were 
measured experimentally for the scattering of p-polarized light from weakly rough, one-dimensional gold surfaces. 
The C^ 15 -*, C^ 2 ', and correlation functions have yet to be observed experimentally. 

The question arises as to whether it possible to determine the relative magnitudes of the different correlation 
functions from a knowledge of the experimental parameters of the surface roughness and its statistical properties. 
This question has been raised earlier in [12,16], but not answered definitively. We therefore address it here for the case 
of a one-dimensional random surface defined by the equation X3 = ((xi), on the basis of the single assumption that 
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the surface profile function C( x i) is a single- valued function of x\ that constitutes a zero- mean, stationary, Gaussian 
random process. 

The outline of this paper is as follows. In Section 2 we introduce the angular intensity correlation function and 
analyze it in terms of the possible statistics of the scattering matrix. In Section 3 we illustrate the conclusions of 
Section 2 for the simple example of the scattering of light from the randomly rough surface of a perfect conductor. 
Finally, in Section 4 we present the conclusions drawn from the results obtained in this work. 



2. THE ANGULAR INTENSITY CORRELATION FUNCTION 



The general angular intensity correlation function C(q, k\q' , k') we study in this work is defined by 

C(q, k\q>, k') = (I(q\k)I(q'\k')) - (I(q\k))(I(q'\k')), (2.1) 

where the angle brackets denote an average over the ensemble of realizations of the surface profile function. The 
intensity I{q\k) entering this expression is defined in terms of the scattering matrix S{q\k) for the scattering of light 
of frequency to from a one-dimensional random surface by 

I(q\k) = ^ (3 \S(q\k)f, (2.2) 

where L\ is the length of the xi-axis covered by the random surface, and the wavenumbers k and q are related to the 
angles of incidence and scattering, 9 and 9 S , measured counterclockwise and clockwise from the normal to the mean 
scattering surface, respectively, by k = (uj/c)sin6 and q = (uj/c)siuO s . In terms of the scattering matrix S(q\k) the 
correlation function C(q, k\q', k') becomes 

C(q,k\q',k') = ^[(S(q\k)S*(q\k)S( q '\k')S*(q'\k')) 

L, x C 

- (S(q\k)S*(q\k))(S(q'\k')S*(q'\k')}]. (2.3) 

Since, due to the stationarity of the surface profile function, (S(q\k)) is diagonal in q and k, (S(q\k)) = 2TrS(q~ k)S(k), 
we introduce the incoherent part of the scattering matrix SS(q\k) — S (q\k) — (S (q\k)) . Then, from the relations between 
averages of the products of random functions and the corresponding cumulant average s [2 1] and omitting all terms 
proportional to 2nS(q — k) and/or 2ir8{q' — fc'), as uninteresting specular effects, Eq. ( |2.3| ) can be rewritten in the 
form 

+ (SS(q\k)SS*(q\k)SS(q'\k')SS*(q'\k')) c }, (2.4) 

where (• • -) c denotes the cumulant average. 

Due to the stationarity of the surface profile function ((xi), (SS(q\k)SS* (q'\k')) is proportional to 2nS(q—k — q' + k'). 
It gives rise to the contribution to C(q, k\q' , k') called (q, k\q' , k'), and describes the memory effect and the 
reciprocal memory effect. Similarly, (SS(q\k)5S(q'\k')) is proportional to 2Tr5(q — k + q' — k') 1 and contributes 



the correlation function C^ 10 \q,k\q',k') to C(q,k\q',k'). The third term on the right hand side of Eq. (2.4) 
(SS(q\k)5S*(q\k)SS(q'\k')SS*(q'\k')) c is proportional to 27r<5(0) = L%, due to the stationarity of the surface profile 
function C( a; i)j an d gives rise to the long-range and infinite-range contributions to C(q,k\q' ,k') given by the sum 
C (1 - 5) (<7, k\q', k')+CW (q, k\q', fc')+C (3) (q, k\q'\k'). Thus, we have separated explicitly the contributions to C{q, k\q' ', k') 
that have been named C^ (q,k\q' ,k') and (q, k\q' , k'). 

What is more, from Eq. ( ]2.4| ) we can easily estimate the relative magnitudes of the different contributions to 
the general correlation function. Indeed, since 2ir5(0) — L%, when the arguments of the delta-functions vanish the 
C^ x \q,k\q' ,k') and C^ 10 \q,k\q',k') correlation functions are independent of the length of the surface L\, because 



they contain [27r<5(0)] 2 . At the same time the remaining term in Eq. (2.4), that yields the sum C^ 1 ^ (q, k\q' , k') 

(q, k\q' , k') + C^(q, k\q'\k'), is inversely proportional to the surface length, due to the lack of a second delta 
function. Therefore, in the limit of a long surface or a large illumination area the long-range and infinite-range 
correlations are small compared to short-range correlation functions, and vanish in the limit of an infinitely long 
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surface. Thus, the experimental observation of the C^' 5 \ C^ 2 \ and C' 3 ' correlation functions requires the use of a 
short segment of random surface and/or the use of a beam of narrow width for the incident field. A detailed discussion 
of the conditions under which they may be observed will therefore be deferred to a separate paper. 

The preceding results are consistent with the usual assumptions and conclusions encountered in conventional speckle 
theory [22,23]. Thus, when the surface profile function is assumed to be a stationary random process, and the random 
surface is assumed to be infinitely long, the scattering matrix S(q\k) becomes the sum of a very large number of 
independent contributions from different points on the surface. On invoking the central limit theorem, it is found 
that S(q\k) obeys complex Gaussian statistics. In this case Eq. (2.4) becomes rigorously 

C(q, k\q', k') = [\{SS(q\k)6S*(q'\k'))\ 2 + \{SS(q\k)6S(q'\k'))\ 2 ] (2.5) 

Liy C 

= C^(q,k\q',k') + C^(q,k\q',k'), (2.6) 

because all cumulant averages of products of more than two Gaussian random processes vanish. The last term on the 
right-hand side of Eq. (2.4) therefore gives the correction to the prediction of the central limit theorem due to the 
finite length of the random surface. 

If it is further assumed, as is done in speckle theory, where the disorder is presumed to be strong, that SS(q\k) 
obeys circular complex Gaussian statistics [22,23], then (SS(q\k)6S(q'\k')) = and the expression for C(q,k\q' , fc') 
simplifies to 

C(q,k\q',k') = ^\(5S(q\k)8S*{q'\k')\ 2 (2.7) 

J-J-y C 

= CU(q,k\q',k>). (2.8) 

This approximation is often called the factorization approximation to C(q, k\q'k') [17]. 

We recall that if the complex random variables F± and F 2 are jointly circular complex Gaussian random variables, 
then the conditions 

(ReF 1 ReF 2 ) = (ImF 1 ImF 2 ), (2.9) 
(ReF 1 ImF 2 ) = -(ImF 1 ReF 2 ), (2.10) 

have to be satisfied. To analyze how the scattering matrix transforms from a complex Gaussian random process into 
a circular complex Gaussian random process we represent the scattering matrix in the form SS(q\k) — SS\(q\k) + 
iSS 2 (q\k). The expressions for the averages of the products of the real and imaginary parts of SS(q\k) can be written 
in terms of (6S{q\k)SS* (q'\k')) and (6S(q\k)SS(q'\k')} 

(<JSi(#)aSi(«'|* / )) = \^ [(SS(q\k)6S*(q'\k')) + (SS(q\k)SS(q'\k'))} (2.11) 

(SS 2 (q\k)SS 2 (q'\k')) = ii?e [(5S(q\k)6S*(q'\k')) - (6S(q\k)SS(q'\k'))} (2.12) 

(SS 1 (q\k)SS 2 (q'\k')) = -\lm [(SS(q\k)6S* (q'\k/)) - (dS(q\k)SS(q'\k'))} (2.13) 

(SS 2 (q\k)SS 1 (q'\k')) = i/m [(SS(q\k)5S*(q'\k')) + (5S(q\k)SS(q'\k'))] . (2.14) 



When q = q' and k — k', the average (5 S (q\k)S S (q\k)) , which is proportional to 2n5(2q — 2k) due to the stationarity 
of the surface profile function, is nonzero only in the specular direction q — k. Therefore, if the surface is infinitely 
long, and if we omit the specular direction, from Eqs. (2.9) - (2.10), and Eqs. (2.11) - (2.14) we see that the scattering 
matrix is a circular complex Gaussian random process. Consequently, apart from the specular direction, the speckle 
contrast p = ^[{(5S{q\k)5S*{q\k)) 2 ) / ' (SS(q\k)SS*(q\k)) 2 } - 1 is unity [22-24]. This result contradicts the well-known 
result of Refs. 23 and 24 that the statistics of the diffuse component of the scattered field is highly non-circular when 
the surface is weakly rough, and only in the limit of very rough surfaces is the circularity of the statistics restored. The 
contradiction stems from the representation of the amplitude of the scattered field as the convolution of a real valued 
amplitude weighting function and the random phase factor in [23,24]. The assumption of a real-valued amplitude 
weighting function, which represents the finite width of the aperture, is identical to the assumption of a finite length 
of the randomly rough surface. As a result, the statistics of the scattering amplitude is nonstationary in [23,24]. In 
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the present work we are interested only in the case where the statistics of the surface profile function, as well as of 
the scattering matrix, is stationary. 

The set of the scattering matrices SS(q\k) is a set of jointly circular complex Gaussian random variables when 
(5S(q\k)6S(q'\k')) vanishes. But when (5S(q\k)SS(q'\k')) vanishes the correlation function C^ 10 ' vanishes, since, 
within a coefficient, C^ w \q,k\q'k') ~ \(SS(q\k)6S(q'\k'))\ 2 . 

Thus, calculations and measurements of the correlation function C(q,k\q',k') yields important information about 
the statistical properties of the amplitude of the scattered field. If the random surface is such that only the and 
£7(10) corre i a tion functions are observed, then S(q\k) obeys complex Gaussian statistics. If the random surface is such 
that only is observed, then S(q\k) obeys circular complex Gaussian statistics. Finally, if the random surface is 
such that C (15) , C* (2) and C (3) are observed in addition to both C^ 1 ** and C^ 10 \ then S(q\k) is not a Gaussian random 
process, but the statistics it obeys in this case are not known at the present time. 

To conclude this section we introduce the the normalized angular intensity correlation functions of interest to us, 
which in terms of SS(q\k) are defined by 

kW y) = \{SS(q\k)SS*(q'\k'W f2 15) 

" [q ' Klq > K) (6S{q\k)SS*{q\k)){SS(q / \k')SS*{q'\k')) , 1 ' 

and 

H(1 o) fo kW k ') = \{6S(q\k)6SW\k>))f 

~ W ' 19 ' > (6S(q\k)5S*( q \k)){6S(q>\k>)SS*( q >\k>))- [ ' °> 
We introduce also the the envelopes Cq 1 ' and Cq 10 "* of the correlation functions and C^ 10 ' , which we define by 

C (1) (q, k\q', k') = 2w5(q -k-q' + k')C^ ] (q, k\q',q'-q + k) (2.17) 

and 

C {W) (q, k\q', k') = 2nS{q - k + q' - k')C^ W) {q, k\q' q 1 + q-k). (2.18) 



3. LIGHT SCATTERING FROM A PERFECTLY CONDUCTING RANDOMLY ROUGH SURFACE IN 
THE FRAMEWORK OF PHASE PERTURBATION THEORY. 

In this Section we study the statistical properties of the scattering matrix for the problem of the scattering of a 
scalar plane wave from a randomly rough infinitely long surface defined by the equation X3 = C( x i)- The region 
2:3 > C( x i) is vacuum, while the region X3 < Ci x i) is a perfectly conducting medium. It is assumed that the Dirichlet 
boundary condition is satisfied on the surface X3 = C(^i)- 

The surface profile function £(xi) is assumed to be a single- valued function of x\ that is differentiable and constitutes 
a zero-mean, stationary, Gaussian random process defined by the properties 

(C(zi)) = 0, (C(^)CK)> = S 2 W(\ Xl - x[\). (3.1) 

In Eqs.( |3.lD the angle brackets denote an average over the ensemble of realizations of ((xi), and S — (C 2 {xi)}^ is the 
rms height of the surface, W(|a;i|) is the surface height autocorrelation function. In numerical examples we will use 
the Gaussian form for W(|a;i|) 

W(\ Xl \) = eM- x l/a 2 ), (3.2) 

where a is the transverse correlation length of the surface roughness. 

A reciprocal phase-perturbation theory for the scattering matrix S(q\k) was constructed in Refs. 25 and 26. The 
term of lowest order in the surface profile function was shown to have the form 

/oo 
d Xl e~ l(q - k)xi e -2V Q o(9)"o(fc)C(*i) _ (3.3) 
-00 

Since 
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(£(#)) = 27T<5(<? - fc) e -2<5 2 «o(9)aoW j 

we can write the expression for SS(q\k) as 

/>oc 

SS(q\k) = / dx ie - l ^^ Xl 



-2i^ a (q)a (k)Q(x l ) _ -25 2 a (q)a (k) 



(3.4) 



(3.5) 



We calculate the averages (SS(q\k)5S*(q'\k')) and (SS(q\k)5S(q'\k')) using the expression ( ftlf ) for the scattering 
matrix. For (5S'(g|A;)5S'*(g''|fc')) we obtain 



(6S(q\k)SS*(q'\k')) 



dx 



e -2i\/ Qo(9)«o(fc)C(^i) _ e -2(5 2 Q!o(9)ao(fe) 



3 2i^/ Q (9') a o(fe')C( a; 'i) _ „-2<5 2 a (g') Q o(fe') 



-25 2 (c to (g)ao(fe)+a («')ao(fc')) 



x 



4<5 2 v /ao(g)Qo(g')"o(fc)ao(fc')W ( | xi—x[ 
2TTS(q -k-q' + fc') e -25 2 (" («)«o(fe)+«o(<?')«o(fc')) 
e 4«5 2 ^a (?)ao(?')"o(fc)"o(fc')M / (l«l) _ 



- 1 



x 



du 



,-i(q'-k')l 



(3.6) 

(3.7) 
(3.8) 



while for (6S(q\k)6S(q'\k r )} we have 



(6S(q\k)5S(q'\k')) = 



dx\ 



dx' 1 e~ i ( q ~ k ^ Xl ~ i ( q ' ~ k ') x 'i 



D -2iy/ a (q)a (k)C(x 1 



-25 a (q)a (k) 



D 2i^a ( 9 ')ao(fc')CK) _ g-2<5 2 a (?' )<*o (*') 



-<5 2 (a (9)ao(fe)+ao(«')ao(fe'))/2 



X 



e -45 2 A /a (g)ao(g') ct o(fe)ao(fe') H/ (l :I: i- a: 'il) _ 1 

27r5(g — jfe + g' — fc')g-2t5 2 (ao(g)ao(fc)+«o(g') Q o(fc')) 
/ r e -4<5V a o(9)«o(9')"o(fe)ao(fe')H'(|«|) _ 



i(q'-k')u 



(3.9) 



(3.10) 



(3.11) 



It is readily seen that in contrast to (SS(q\k)SS* (q'\k')} the average (SS(q\k)SS(q'\k')} vanishes with incr easin g 
roughness parameters S and a, due to the negative exponential under the integral sign in the last line of Eq. fl3.1l| ). 
Plots of the normalized correlation functions SW(<7, A;|g', A:') and S( 10 )(g, k') as functions of <5 for different values 
of a are presented in Fig. 1 (a), while plots of the envelopes of the correlation functions and C' 10 ) as functions 
of 8 for different values of a are presented in Fig. 1(b), for fixed values of q, k and q' , while k' is determined by the 
constraint of the corresponding 6— function. When calculating the results presented in Figs. 1 (a) and (b) the value 
of q' was chosen to produce the same values of and C^ 10 ^ in the limit of a weakly rough surface. From the plots 
presented in Fig. 1(a) we see that ^ 10 \q, k\q' , k') vanishes even for quite moderately weakly rough surfaces for which 
E^(q, k\q', k') is still about unity. We note that C^ 1 ' also decreases with increasing 5 (Fig. 1 (b)). Using Eqs.(2.11) 
- (2.14), Q and fyl) we obtain the expressions for {{SS^k)) 2 ) , {{5S 2 (q\k)f) and (SS 1 (q\k)5S 2 {q\k)) 



((SSxiqlk)) 2 ) = e -45 2 ao(9)oo(fc) 

+ -7T% - k) 



2 



du cos 



(q - k)u 



S 2 a (q)ao(k)W(\u\) 



>s(q — k)u C 



e -8 2 ao (q) aa (k)W(\u\) _ 1 



(3.12) 



and 



5 



1.20 



~4 0.f 



0.40 



0.20 



0.00 



0.00 



a=300nm 
a=500nm 
a=800nm 



0.20 



0.40 



0.60 



m. 



0.10 



x 

I 

cr 

I 

q 

"a 
x. 
u 



0.05 



0.03 



; a 


(b) 


/ \ 

- ,,\ 


a=300nm 




— - a=500nm 




a=800nm 



0.00 



0.20 

ill 



0.40 




0.08 



a=300nm 
a=500nm 
a=800nm 



2 0-05 
i 

u 

I 

u 

I 0.03 





(d) 




-a=300nm 




a=500nm 




a=800nm 


// ■ x25 \ \ 





0.00 



0.20 

ill 



0.40 



FIG. 1. The normalized correlation functions (a) and H' 10 - 1 (c), and the envelopes Cq 10 ' (b) and Cg 10 " 1 (d) as functions 
of 8/ A for values of the transverse correlation length a — 300 nm, 500 nm, and 800 nm. The incident light was s— polarized 
and of wavelength 632. 8nm. The scattering medium was a randomly rough perfect conductor. Furthermore do = 30°, 9 S = 0°, 
and 6' s — 0° . In Fig. la the results for the different correlation lengths considered could not be distinguished. 



{{5S 2 {q\k)f)=e 



-4S 2 a (q)a (k) 



h, 

2 



/oo 
ducos(q - k)u / e <5 2 «o(g)ao(fc)W(M) _ A 
-oo 

1 f'°° / \ 

-nS{q -k) J cos(q - k)u ^ 2 ^(^o(k)W(\u\) _ ^ 



while 



(SS 1 (k\k)6S 2 (k\k)) = 0. 



(3.13) 



(3.14) 



In Fig. 2 we present plots of the ratio ((5S2(k\k)) 2 } / ((SSi(k\k)) 2 ) as a function of the rms height of the surface 
roughness 5. Since this ratio is calculated for the specular direction q = k, it is independent of the transverse 
correlation length a. From the plot presented it is easily seen that for large values of the rms height the incoherent 
part of the scattering matrix, SS(q\k), becomes a circular complex Gaussian variable, even in the specular direction. 
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FIG. 2. The ratio ((SS 2 (k\k)) 2 ) / ((SSi(k\k)f) as a function of 6/X. 



4. LIGHT SCATTERING FROM A RANDOMLY ROUGH PENETRABLE SURFACE 



The results of the preceding Section enable us to make several conclusions when studying the scattering of light 
from a randomly rough surface of a penetrable medium. For simplicity we consider here the scattering of s— polarized 
light from randomly rough surface of a medium characterized by a dielectric function e(u>). As is well known (see, e.g. 
[27-29]) if the surface profile function is such that the conditions for the applicability of the Rayleigh hypotesis are 
satisfied the scattering amplitude R{q\k) obeys the reduced Rayleigh equation. Rewritten in terms of the scattering 
matrix S(q\k) it has the form 



S(q\k) = 2n5(q - k)R (k) + N(q\k) + J°° ^-M(q\p)S(p\k), 



where for the case of the scattering of s— polarized light, 

ao(k) — a(k) 



R (k) 



cto{k) + a(k) ' 



(4.1) 



(4.2) 



a (k) 



k 2 , a(k) = \ e{uj) 



fc 2 , 



(4.3) 



N(q\k) = -{e-l) 



LU 2 /C 2 



<*o(q) J(a{p) + a (k)\p - fc) 



a (q) + a{q) V a {k) a{p) + a (k) 



(4.4) 



and 



M(#) = -(e-l)) 



la (q) J(a(p) - a (k)\p- k) 



a (q) + a(q) V a (p) a(p) - a (k) 



J{j\Q) = f dx ie - lQxi ( e ~ l ^ {xi) - l) 

J — oo 



(4.5) 



(4.6) 
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FIG. 3. The envelopes of the C 1 ' 1 ' (a) and C' 10 ' (b) correlation functions as functions of 6' s for 60 = 30° and 9 S = 0° , while 
0' is determined by the constraints of the corresponding S— functions, for the scattering of s— polarized light from a randomly 
rough silver surface with a — 500nm and S = 20 nm (solid lines), S = 50 nm (dashed lines), and 5 — 100 nm (dotted lines). 



We can write the solution of Eq. (4.1) formally as 



S{q\k)=Ro(k)2irS(q-k)+F(q\k)+ I ^M(q\p)F(p\k) 



2tt 

+ /I ^ M{qlp) /I (PbWl*) + ■■■, (4.7) 

where 

F(q\k) = N(q\k) + M(q\k)R (k), (4.8) 

and we keep all terms in the infinite iterative series. Both N(q\k) and M(q\p) contain the surface disorder only in the 
functions J(p/\Q). Therefore, having in hand the recipe for calculating the average of the product of any number of 
functions J{^\Q), we can calculate, in principle, both (5S(q\k)5S(q'\k')) and (SS(q\k)SS*(q'\k')). The basics of such 
calculations were described in Ref. [30]. 

To calculate the averages (SS(q\k)8S(q'\k')) and (8S(q\k)SS*(q'\k')) we multiply the series (4.7) for S(q\k) by the 
corresponding series for S(q'\k'), and average the product term-by-term. From the result we subtract the product 
(S (q\k)) (S (q' \ k')} . In a similar fashion we calculate the average (5 S (q\k)S S* (q' \ k')} by multiplying the series (4.7) for 
S(q\k) by the complex conjugate of the corresponding series for S(q'\k'), averaging the product term-by-term, and 
subtracting the product (S(q\k))(S*(q'\k')) from the result. In the product(<5 S(q\k)SS*(q'\k')) the contribution of n th 
order in the functions J(7|Q) and J*(7|Q) contains n — 1 terms of the form 

n— 1 ( J m n—m \ J m \ ln — vn \ ^ 

e ( n j ^\Qr) n j*Ws) ) ( n ) i n j*«\q' s ) ) ■ ^ 

m=l I \r=l s=l / \r=l I \ s=l I ) 

To obtain a nonzero contribution, for each value of m at least one J("f r \Qr) must be contracted with at least one 
J*(Y S \Q' S )- Therefore, each term in this sum contains at least one factor with a postive exponential of the form 
exp{<5 2 77'lF(|u|)} — 1. In contrast, when calculating (SS(q\k)SS(q'\k')) the contribution of the n th order in the 
functions J(j\Q) contains the sum 

n— 1 ( J m n—m \ I m \ /n — m \ \ 

E n j MQ r ) n j «\q's) ) ( n ) i n j «\q's) ) ■ ^ 

ra=\ \. \r=\ s=l / \r=l / \ s=l / ) 
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FIG. 4. The same as in Fig. 3, but for a = 3.85/im and 5 — 1.278/im (solid lines) and 8 = 0.1278/im (dashed lines). 



In this case, to obtain a nonzero contribution, for each value of m at least one J("/ r \Qs) must be contracted 
with at least one J(j' s \Q' s ). Therefore, each term in this sum contains only negative exponentials of the form 
exp{— <5 2 77'W / (|u|)} — 1. Owing to this lack of the positive exponential, (SS(q\k)SS(q'\k')) vanishes when the roughness 
parameters increase. 

In Fig. 3 we present plots of the envelopes Cq 1 ^ and Cg 10 ^ of the corr elation functions (Fig. 3(a)) and (Fig. 
3(b)) as functions of 9' s for fixed values of 9q and 9 S , while 9' Q is determined by the constraints of the corresponding 
5— functions. The calculations were carried out for the scattering of s— polarized light, of 612. 7nm wavelength, from 
a weakly rough random surface of a silver characterized by the complex dielectric constant e = —17.2 + i0.479 for 
different values of the roughness parameters 5 and a. In calculating the results presented in Fig. 3 we kept all terms 
in the infinite iterative series Eq. (F7) which would give the contributions to the averages we calculate through terms 
of 0(S 8 ) if they were to be expanded in powers of the small parameter (lo/c)S. 

In Fig. 4 we present rigorous numerical simulation calculation results |3l| for the envelopes of the correlation 
functions C' 1 ' (Fig. 4a) and C^ 10 ' (Fig. 4b). The surface parameters used here were the same as those used in 
obtaining Fig. 3 except that the roughness now was 8 = 1.278 /im (solid lines) and 5 — 0.1278 fj,m (dashed lines). 
It should be pointed out that for the scattering of s-polarized light from a weakly rough random metal surface 
there should be no memory- or reciprocal memory-effect present in . This is indeed confirmed by our numerical 
calculations where the Cq for S — 0.1278 \xm (Fig. 4a, dashed line) is a smooth function of its argument, as well 
as by the results presented in Fig. 3a. In particular, there are no peaks at angles 9 — 0° and 30°, which are the 
positions of the memory- and reciprocal memory-effects. As the roughness is increased to S — 1.278 /im one sees from 
Fig. 4a (solid line) that the overall amplitude of the envelope Cq 1 ' is increased and, more important, that two peaks 
have developed at the aforementioned angles. These peaks are due, in the large roughness limit, to volume waves 
scattered multiply at the rough surface. In Fig. 4b the corresponding results for the Cq 10 ^ -envelopes are presented. It 
is observed that in the low roughness limit this envelope is structureless, and that and C$ 10 ^ are roughly of the 
same order of magnitude. However, as 6 is increased, the scattering matrix S(q\k) starts to obey circular complex 
Gaussian statistics, and thus as discussed earlier, the envelope Cq 10 ^ should in principle vanish. From our numerical 



results for 6 = 1.278 [im (solid line) we indeed see that Cq 10 "* is much smaller then the corresponding Cq 1 ' shown in 
Fig. 4a. In fact Cq 10 ^ is just noise, consistent with this function vanishing in the large roughness limit. 



.(i) 
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5. CONCLUSIONS 



In this paper we calculated the angular intensity correlation functions by means of an approach that explicitly 
separates out different contributions to it. We have shown that calculations and measurements of the correlation 
function C{q, k\q' , k') yields important information about the statistical properties of the amplitude of the scattered 
field. In particular, we have shown that the short-range correlation function C'^ 10 ^ is, in a sense, a measure of the 
noncircularity of the complex Gaussian statistics of the scattering matrix. Thus, if the random surface is such that 
only the and C^ 10 ) correlation functions are observed, then S(q\k) obeys complex Gaussian statistics. If the 
random surface is such that only is observed, then S(q\k) obeys circular complex Gaussian statistics. Finally, if 
the random surface is such that C^' 5 \ and C^ 3 ' are observed in addition to both and C^ 10 \ then S(q\k) is 
not a Gaussian random process. 
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